Breakdown of the classical description of a local system 
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We provide a straightforward demonstration of a fundamental difference between classical and 
quantum mechanics for a single local system; namely the absence of a joint probability distribution of 
the position x and momentum p. Elaborating on a recently reported criterion by Bednorz and Belzig 
[Phys. Rev. A 83, 52113] we derive a simple criterion that must be fulfilled for any joint probability 
distribution in classical physics. We demonstrate the violation of this criterion using homodyne 
measurement of a single photon state, thus proving a straightforward signature of the breakdown of 
a classical description of the underlying state. Most importantly, the criterion used does not rely on 
quantum mechanics and can thus be used to demonstrate non-classicality of systems not immediately 
apparent to exhibit quantum behavior. The criterion is directly applicable any system described by 
the continuous canonical variables x and p, such as a mechanical or an electrical oscillator and a 
collective spin of a large ensemble. 
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The conceptual differences between classical and quan- 
tum physics have intrigued and sometimes bewildered the 
physics community since the early days of quantum me- 
chanics. This has led to a search for indisputable man- 
ifestations of the quantum world through observations 
of non-classical behavior in experiments. A field of par- 
ticular curiosity is that of identifying the quantum to 
classical cross-over for ever larger systems, thereby even- 
tually identifying non-classical effects in macroscopic sys- 
tems. Recently this has led to the observation of, e.g., 
macroscopic entangled atomic ensembles [l|, Q , interfer- 
ence of large molecules [3|] and experiments pushing to- 
ward observing non-classical effects in mechanical oscil- 
lators [J^- In parallel to this fundamental interest, non- 
classicality is of central importance to quantum informa- 
tion processing, the essence of which is to advance com- 
putation beyond what is classically possible However, 
in some instances quantum effects are claimed by demon- 
strating consistency with an appropriate quantum model. 
Yet any rigorous demonstration of genuine quantum be- 
havior must exclude the possibility of classical explana- 
tions. The importance of this is exemplified in Ref. Q, 
where a pair of coupled classical oscillators is shown to ex- 
hibit signatures easily mistaken for those of entanglement 
expected from a quantum model. Thus, a definite conclu- 
sion on the quantum nature of a system can only result 
from the breakdown of the classical description and not 
from verified agreement with quantum mechanics. This 
approach is most rigorously demonstrated by the Bell- 
inequalities, where the underlying model of the system is 
stripped of any physics and is reduced to the very basic 
assumptions of locality and realism, resulting in an indis- 
putable non-classicality criterion. The Bell-inequalities 
cannot however, by their very nature, be investigated by 



data obtained from a single system. 

In this letter we provide a conceptually simple demon- 
stration of one of the key discrepancies between classical 
and quantum mechanics, valid for systems of a single 
degree of freedom: classical systems can always be de- 
scribed by a joint probability distribution for x and p, 
the two canonically conjugated coordinates of a system, 
whereas such a description does not apply in quantum 
mechanics due to the Heisenberg uncertainty principle. 
This discrepancy is most evident when the phase space 
description of the state of a system is examined. Classi- 
cally, the phase space distribution W{xi,pi) is the joint 
probability of finding the system in an infinitesimal area 
around x = Xi,p = pi, and hence it obeys all the require- 
ments of a probability distribution including being a non- 
negative function. As mentioned, in the case of a quan- 
tum phase space formulation, introduced by Wigner 
the Heisenberg uncertainty renders this definition mean- 
ingless, as a joint probability distribution for x and p 
does not exist. The phase space distribution is only de- 
fined through the single coordinate (marginal) distribu- 
tions, projected from the distribution function jXl] and 
this relaxation of constraints allows for negative values 
of the function in areas smaller than h. This negativity 
is not directly observable due to the vacuum fluctuations 
preventing simultaneous measurement of x and p. How- 
ever, one can still infer the phase space distribution from 
measurements of only a single observable at a time and 
detect such negativities, thereby illuminating the failure 
of classical theory. 

The usage of these negativities as markers of non- 
classicality has been discussed and demonstrated in sev- 
eral quantum optics systems (see, e.g., (lli - [l5| ). using 
tomographic techniques. Often such methods search 
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for the quantum state most compatible with the ex- 
perimental data using statistical inference or variational 
techniques [l6. 17 1 and thus inherently rely on quantum 
mechanics. These methods are therefore not applicable 
for demonstrating the absence of a classical description. 
Alternatively, given measurements of all the coordinate 
distributions, the underlying state can be uniquely de- 
termined, and the phase space distribution fully calcu- 
lated using the inverse Radon transformation [18[ with- 
out relying on quantum mechanics. Though such meth- 
ods have been used in quantum optics for demonstrating 
various states, the mathematical transformation involved 
is highly complicated. Furthermore the numerical stabil- 
ity of the inverse transformation is problematic, leading 
to numerical uncertainty at high frequencies, and some- 
times results in unphysical states These limitations 
are a drawback for using tomographic techniques for val- 
idating the breakdown of a classical description, and the 
application of these methods is usually cumbersome. 

Our simple, unambiguous demonstration of the ab- 
sence of a classical probability distribution is based on 
recent theoretical work by Bednorz and Belzig \2d\ that 
verifies the negativity of the Wigner function based on 
moments. As discussed in detail below, their results lead 
to a hierarchy of inequalities, such that violation of any 
one inequality indicates negativity of the Wigner func- 
tion. Full tomographic reconstruction with the associ- 
ated numerical complexities is thereby avoided. We ex- 
tend this approach such that it can be applied to quadra- 
ture measurement of a single photon state, and use the 
experimental data from the heralded single photon gen- 
eration to directly disprove the existence of a joint prob- 
ability of the position and momentum for this system. 

We start by re-iterating the key results of Bednorz and 
Belzig, through a reformulation that relies only on clas- 
sical mechanics. The phase space of a system with a 
single degree of freedom is fully characterized by a two- 
dimensional phase space distribution W{x,p). That is, 
given the phase space distribution, the ensemble averaged 
result of any measurable quantity A can be obtained by 



(A) 



dxdp W{x, p) A{x, p) 



(1) 



where A(x^p) is the decomposition of the quantity A in 
terms of the generalized coordinate x and its canonically 
conjugated momentum p. 

To disprove the existence of a classical probability dis- 
tribution we examine the ensemble average of a non- 
negative test function ^{x,p) over a classically explain- 
able system, which must have a proper distribution func- 
tion that results in the ensemble average of 5 be non- 
negative: 

(5) - j dxdpW{x,p)S{x,p) > 0. (2) 
Violating this condition is a direct proof of the absence of 



a joint probability distribution. The condition can, how- 
ever, be violated in quantum mechanics, where W{x,p) 
is the Wigner function that can contain negative values. 
The objective therefore is to optimize a test function such 
that it will be dominant at the possible negative areas of 
the distribution function. For a rotationally invariant 
phase space both the phase space distribution and the 
test function can be described solely by the phase space 
radius r, defined by = x^ +p'^. For reasons to be- 
come clear later, we choose a specific form for the test 
function writing it as a square of an A'^th order, even 
polynomial M with real coefficients {Ci}jv; 



id) = {M') = 




(3) 



Minimizing the above expression for a given order A^ is 
done by straight-forward linear optimization of the coef- 
ficients {Ci}^- 



N/2 



(4) 



for all j = 1,2, N/2. Notice that the linearity of the 
problem ensures that the obtained minimum of (J?) is 
global and therefore the most optimal indicator of a pos- 
sible violation of Eq. ^ for a given polynomial order 
N. It is important to emphasize that this is only a suffi- 
cient criterion for non-classicality, and an optimized pos- 
itive average for a chosen N does not ensure a classical 
probability distribution, since the negativity may only 
be exhibited by the inclusion of higher order terms in 
M. However, it is clear that increasing the polynomial 
order A^ cannot increase the minimized value of {^) , and 
we conjecture that the limit oi N ^ oo will exhibit any 
negativity of the Wigner function, as the polynomial can 
represent an arbitrarily (analytical) sharp peaked func- 
tion 5 focused at the negativity. Assuming the existence 
of all moments (e.g., due to an exponentially decaying 
tail of the phase space distribution at large r) , this then 
becomes a necessary criterion for the negativity of the 
distribution function. We also note here, that similar 
polynomial expansion has been discussed 21, 2^ in the 
context of the P-function distribution. The P-function 
is, however, only defined within the framework of quan- 
tum mechanics, and hence cannot be used to prove the 
absence of a classical description. 

We assume that, as is the case for many systems, the 
system in question can only be experimentally accessed 
by measuring one of the canonically conjugated variables 
(e.g., X or p) at a time. Since we are restricted to single 
coordinate measurement at a time, neither the intensity 
nor the phase space distribution function is directly ac- 
cessible. For this method to be applicable to such ex- 
perimental data, the functional {^) must be expressed in 
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terms of the moments of the projected coordinates (Q"), 
where 

Qa = cosa X + sinap, Pa = cosax — sinap (5) 

is a measureable rotated coordinate. To do this we use 
the identity 

2N 2Ar 

(a,2+p2)^ ^ ^2^^ (cos(l^) x + p) ,(6) 

m— 1 

where 

This is where quantum and classical approaches diverge. 
While classically Eq. ^ represents a measurable phys- 
ical quantity, it is missing the key vacuum uncertainty, 
allowing for the breakdown of the classical description. 

It is interesting to note the implication of identity ^ . 
For the 2mth moment of the radial distribution to be 
known, we need 2m 'cuts' in phase space; i.e., different 
coordinate measurements at equally distributed angles. 
Regardless of any assumption about the underlying state, 
the average of Eq. ^ directly gives 

2N 

(r2^)=A2A.5]((Q_/2A.)'^). (8) 
m— 1 

In the special case of a symmetric distribution function 
these moments are all identical, and Eq. ([8|) reduces to 

{r'^)=^^l^y\'^{x'^). (9) 

The radial moments can thus be indirectly calculated 
from the quadrature measurements. Substituting these 
radial moments into Eq. ^ using the functional form of 
d{x,p) given by Eq. ([3]), we get, for a given set of mea- 
sured moments {{x'^'')}k, a necessary condition for classi- 
cality of the underlying state. If Eq. ([2|) is violated by the 
solutions of Eq. (jU, the underlying state cannot be ex- 
plained by a proper phase space probability distribution, 
and one cannot assign a joint probability distribution to 
X and p. 

To demonstrate the absence of a joint probability dis- 
tribution we are going to consider the phase space de- 
scription of a single photon state. In phase space this can 
be described by the first excited state of a harmonic oscil- 
lator, which is rotationally invariant and contain negative 
parts in the Wigner functions. Fig. [1] shows the optimal 
functional forms obtained for the this state for low poly- 
nomial orders. As higher order terms are included, the 
optimized test function is increasingly probing the nega- 
tive part of W, yielding a negative expectation value. We 
note that negative expectation values appear only from 
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FIG. 1: Profiles of the test function ^ minimizing the ex- 
pectation value (5) for the first excited state of a quantum 
harmonic oscillator, as a function of the phase space radius, 
for different orders A'^ (see text) plotted against the profile 
of the corresponding Wigner function. As the order of the 
polynomial increases, the function becomes centered around 
the negativity, decreasing elsewhere. In this case, negative ex- 
pectation values are obtained starting at A'^ = 4. Inset shows 
the polynomial order required to observe negative expectation 
values, as a function of the single-photon fractional content 
in a mixture with vacuum. As the fraction of vacuum is in- 
creased, the state approaches a classically describable state 
and higher moments are needed to observe the negativity. 



the fourth order onward. This is because the peak of 
the test function at the position of the negativity must 
be narrower than Heisenberg's uncertainty in order not 
to smear the negativity; this is in full agreement with 
Ref. [13. 

The experimental demonstration is achieved with sin- 
gle photons generated by an heralded cavity-enhanced 
non-degenerate parametric down-conversion. The equiv- 
alence between a single mode electromagnetic field and 
an harmonic oscillator allows us to describe the EM field 
by a phase space of a single degree of freedom. The down- 
conversion process produces two photons, and as one is 
detected as a trigger, the result is a single photon state 
where the losses introduce a statistically mixed compo- 
nent of vacuum. The projection measurements (quadra- 
tures) are obtained by measuring the statistics of the 
noise, using an optical homodyne detection scheme. In 
this scheme, the weak investigated optical field is over- 
lapped with a strong laser pulse on a beam splitter, and 
the interference of the two fields is detected and sub- 
tracted. The phase of the strong laser field determines 
the angle a (Eq. ([5])) of the measured coordinate. Mea- 
surements were taken without fixing the phase of the lo- 
cal oscillator, thus smearing the resulting distribution. 
This enables us to treat the results as rotationally invari- 
ant even if non-invariant features existed prior to smear- 
ing. Such measurements will generate a rotationally in- 
variant reconstructed state for any underlying state, but 
this does not necessarily average out negativities in the 
Wigner function (see, e.g., [l^). For details of the ex- 
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FIG. 2: Expectation value for the square of a polynomial 
relative to its standard deviation, as a function of the polyno- 
mial's order for the experimental data. Negativity by almost 
20 standard deviations disproves the existence of a joint prob- 
ability distribution for x and p. The inset shows a histogram 
of the raw measured quadrature data (arbitrary units) . 



perimental setup and the charactarization of the resulting 
single photon see Ref. j24l] . 

The data set contained 180,000 measured quadratures. 
We have here revised the optimization of the functional 
to also account for statistical uncertainties inherent to a 
limited data set. This is done by optimizing 
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where cry = \/ (5^} — (i?)^ is the standard deviation of 
^. The results are shown in Fig. [21 The fact that the 
expectation value for our test function is negative with 
certainty of almost twenty standard deviations, clearly 
demonstrates that the measured state in this experiment 
cannot be explained by classical theory, unambiguously 
negating the possibility of existence of a joint probabil- 
ity distribution for x and p. The appearance of negative 
values from the twelfth order polynomials and higher in- 
dicate the quantum mechanical description of this state 
in terms of a Wigner function includes negative valued ar- 
eas. We note that the minimized function from Eq. ([3|) is 
monotonically decreasing for increasing order N, and the 
onset of negativity at a certain order therefore means that 
all higher orders will also be negative. This suggests a se- 
quential authentication procedure for an unknown state. 
As mentioned above, for a pure single photon state, neg- 
ative expectation values are observable from the 4-th or- 
der polynomial onwards. The twelfth order polynomial 
required here is due to the vacuum component of the 
field, requiring higher orders of the polynomial as shown 
in the inset of Fig. [2land is in agreement with the re- 
sults obtained in Ref. 2J] reporting 62% fraction of single 
photon in the resulting mixed state. 

In conclusion, we have experimentally demonstrated 
the non-existence of a joint probability distribution of 
two canonical variables. This is done by violation of an 



inequality derived without the assumptions of quantum 
mechanics, thus allowing for it as proof of the absence of a 
classical description in systems not immediately evident 
to display quantum behavior. The procedure used here 
can thus provide a simple, practical tool for demonstrat- 
ing the non-classicality of a state based on quadrature 
measurements, where the existence of a classical joint 
distribution of two conjugated variables can be negated. 
In this way, this procedure is closely linked to other crite- 
ria 25l-l27| demonstrating contextuality of measurements. 



and thus disproving the classical local hidden variable 



view. Unlike Ref. |25l - l27l |. which are applicable to dis- 
crete variables, the method demonstrated here applies 
for continuous variables such as position and momen- 
tum, collective spin operators and quadrature phase 
operators. This makes it useful to systems containing 
many particles, where criteria based on counting parti- 
cles are not easily implemented and interpreted. This 
method complements the full tomographic reconstruc- 
tion techniques in that it is simpler and avoids numerical 
complexities of inverse transformations. These kinds of 
conceptual proofs, when extended to different detection 
schemes, can shed more light on the quantum to classical 
correspondence, especially where the control of claimed 
macroscopic quantum states is in question. 
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